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1. Introduction 



One of the outstanding problems in the theory of gravitation (and more generally in the 
quest for a unified theory of all interactions) is the singularity problem. According to the 
Penrose-Hawking theorems 1 ), general relativity (GR) manifolds are, in general, geodesically 
incomplete, which is a sign that singularities in space-time occur. 

Singularities are undesirable for a theory which claims to be complete, since their ex- 
istence implies that space-time cannot be continued past them. The space-time structure 
becomes unpredictable already at the classical level. 

Two important examples of singularities in GR are the initial and final singularities in 
a closed Universe and the singularity in the center of the black hole. In the former case, the 
singularity implies we cannot answer the question what will happen after the "big crunch" , 
or (in the case of an expanding Universe) what was before the "big bang" . 

The presence of singularities is an indication that GR is an incomplete theory. Wheeler 
even talks about a "crisis in physics" 2 -*. It is a widespread opinion that either quantum 
gravity or a more fundamental theory such as string theory will provide a cure for the 
"sickness" of GR. However, quantum gravity does not yet exist as a self-consistent non- 
perturbative theory. Neither does string theory exist as a unique theory capable of addressing 
the singularity problem of gravity in a definitive way, although interesting string-specific 
ideas have recently been put forwards 3 ). 

Because of the absence of a completely developed fundamental theory on the basis of 
which we could address the singularity problem, we will use a rather different approach. 
Any fundamental theory will, in the region of low curvature, give an effective action for a 
four-dimensional space-time metric g^y which to lowest order must agree with the Einstein 
action. We will try to construct (guess) an effective action for g^ v which solves the singularity 
problem and which in the low curvature limit reduces to the Einstein action. It is possible 
that in such a manner we will be able to discover important features of the future fundamental 
theory. We might also gain information which will help in finding this fundamental theory. 

Before discussing the ideas behind our construction of the effective action for gravity, 
we return to the Penrose-Hawking theorems 1 ). They do not give us any detailed informa- 
tion about the nature of the singularity. However, in the two examples discussed above - 
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collapsing Universe and black hole - we know that at the singularity some of the physically 
measurable curvature invariants like R, R^ v R^ v , C 2 = C a p 7 s c a ^ s diverge (here, R is 
the Ricci scalar, R^ v the Ricci tensor, and C a ^ the Weyl tensor). It is reasonable to as- 
sume that the divergence of some curvature invariants at the singularity is a fairly general 
phenomenon.. In fact, for singularities reached on timelike curves in a globally hyperbolic 
space-time it can be proved 4 ) that the Riemann tensor becomes infinite. Hence, as a first 
step we will find a mechanism to bound all the curvature invariants. 

Limitation principles play a very important role in physics. Special relativity includes as 
one of its fundamentals the principle that no particle velocity can exceed the speed of light. 
The cornerstone of quantum mechanics is the uncertainty principle which states that the 
second fundamental constant, Planck's constant h, gives the minimal phase space volume a 
particle can be localized in. The third fundamental constant, Newton's gravitational constant 
G, has not yet been used in any limitation principle. 

1 /9 

Thus it is natural to assume that there exists a fundamental length £ p £ ~ (GTic^) ~ 
1CP 33 cm in nature (determined by G) such that there is no curvature corresponding to scales 
i < £ p £. There are strong indications that this will in fact arise in quantum gravity 5 -* or 
string theory 3 ^ . From the existence of a fundamental length it follows by simple dimensional 
considerations that all curvature invariants are limited: 

\R\ < L-} , \R^RH < t£ , \C aPl5 C a ^ s \ <i~ t \... (1.1) 

To realize the idea of a fundamental length it is necessary to construct a theory in which all 
curvature invariants are bounded. Since there are an infinite number of curvature invariants 
and since bounds on low order invariants do not necessarily imply bounds on higher order 
invariants, it is a rather formidable task to construct such a theory. 

Fortunately we can simplify the problem drastically by making use of the "limiting 
curvature hypothesis" (LCH) construction 6 ), according to which one looks for a theory in 
which: 

i) A finite number of invariants are bounded by an explicit construction (e.g. \R\ < £ p 2 
and \R^R^\ <e~l). 

ii) When these invariants take on their limiting values, any solution of the field equations 
reduces to a definite nonsingular solution (e.g. de Sitter space). 
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In this case it follows automatically that all curvature invariants are limited. 

The LCH incorporates in a natural manner Penrose's hypothesis 7 ) that the Weyl tensor 
C should vanish at the beginning (and end) of the Universe. This follows since by the LCH 
the Universe near the big bang and big crunch is de Sitter and that C = for a de Sitter 
Universe. 

A theory in which the LCH is realized has some attractive features, both for cosmology 6 ) 
and for black holes 8 ). In cosmology, the present homogeneous expanding Universe would 
have started out with a de Sitter phase. In this case, we would have some (maybe unusual) 
realization of the oscillating Universe scenario. Entropy considerations tell us that only 
for a perfectly homogeneous and isotropic Universe could we have perfect periodicity. In 
general, we must have a nontrivial realization. Including inhomogeneities, we might obtain a 
multi-Universe model in which one collapsing Universe splits into several de Sitter bounces. 

For black holes, the LCH gives the attractive picture that inside of the horizon instead of 
a singularity at the center we would have a piece of a de Sitter Universe which could be the 
source of other Friedmann (baby) Universes (see Fig. 1). In this case, the difficult question 9 ) 
concerning information loss when matter falls into a black hole have a natural answer: the 
information which is lost to the observer external to the Schwarzschild horizon is stored in 
the baby Universe. In addition, using this picture provides a good starting point to attack 
the issue of the final stage of an evaporating black hole, a problem which has recently been 
of high interest in the context of two-dimensional quantum gravity 10 ) . 

In this paper we construct an effective action for gravity in which all homogeneous and 
isotropic solutions are nonsingular and at high curvature approach de Sitter space. (A brief 
summary of our work was published in Ref. 11.) In order to implement the LCH we procede 
in analogy to a technique by which point particle velocities can be limited, thus achieving 
the transition between Newtonian mechanics and point particle motion in special relativity 
(SR) (see also Ref. 12). An extension of our construction to inhomogeneous cosmologies and 
to black hole metrics will be presented separately. 13 ) 

In the following section we present the general theory of how to implement the LCH. We 
obtain a fairly general effective action for gravity as a higher derivative modification of the 
Einstein action, specialize to the case of an isotropic, homogeneous Universe and derive the 
resulting equations of motion. 
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In Section 3 we analyze a simple model which yields a nonsingular Universe without 
limiting curvature. We discuss the effects of including spatial curvature (i.e. k ^ 0) and 
hydrodynamical matter. The analysis of the more complicated model with limiting curvature 
is given in Section 4. Section 5 contains conclusions and further discussion. 

2. Theory 

In order to realize the LCH and hence to avoid singularities it is necessary to abandon 
at least one of the key assumptions on which the Penrose-Hawking theorems are based. The 
two most important assumptions are: 

i) The energy dominance condition, a simplified version of which appropriate for cosmol- 
ogy is e > and e + 3p > 0, where e and p are matter energy density and pressure 
respectively. 

ii) The Einstein equations are universally true. 

There is no reason to believe that these assumptions will be valid at very high energies and 
curvatures. First of all, already in matter theories routinely studied by particle physicists, 
the energy dominance condition is not always true. For example, the effective equation of 
state for a homogeneous, slowly varying scalar field configuration with potential energy is 
p ~ — e, thus violating the energy dominance condition. This matter evolution scenario is in 
fact the basis for the inflationary Universe 14 ). 

Note, however, that inflationary Universe models do not cure the problem of the final 
singularity. There may be nonsingular solutions for a collapsing Universe filled with scalar 
field matter, but they are of measure zero. Rather, in this case the typical solutions have an 
effective equation of state p = e (the kinetic term for the scalar field dominates), not p = — e, 
and hence have a final singularity. Our goal is to construct a theory in which all solutions 
are nonsingular. 

Concerning the second key assumption of the Penrose-Hawking theorems, it is well 
known that the Einstein theory can only be the effective theory of gravity at low curvatures. 
Perturbative quantum gravity calculations 15 -* , vacuum polarization effects of quantum matter 
fields in an external gravitational background 16 -* , and also considerations based on string 
theory 17 ) all show that the effective equations for the gravitational field should be modified 
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at higher curvatures. In a perturbative analysis, the modifications take the form of higher 
derivative terms which are usually important only at very high (Planck) curvatures. Hence, 
provided the effective action approach is valid at all at high curvatures, this effective action 
will certainly not be of pure Einstein form. 

To summarize, there are two ways to modify the theory at high curvatures in order to 
avoid singularities: 

i) Modify the matter action by including terms which violate the energy dominance 
condition. 

ii) Modify the gravitational field equations. 

The first approach was explored in Ref. 18. However, the weakness of this approach 
is the absence of a good physical motivation for the modification. In addition, it seems 
impossible to avoid singularities associated with purely gravitational modes which do not 
couple to matter. 

The second approach is much better motivated since higher derivative correction terms 
to the Einstein action are predicted by many theories 15-17 ). Hence, our starting point will 
be to look for an effective action for gravity of the form 

S g = - -J^ f F (R , R^BT , C a ^ s C a ^ S , . . .) V=9d 4 x + nonlocal terms , (2.1) 

where the dots denote the dependence of F on other curvature invariants. At low curvatures, 
the leading term in F is simply R. 

The action (2.1) can be viewed as the effective action of some fundamental theory such 
as quantum gravity or string theory. In these theories we are at present unable to calculate 
the nonperturbative effective action. Hence, as mentioned in the introduction, our approach 
will be to construct (guess) an effective action of the form (2.1) to obtain a theory in which 
all solutions are nonsingular. 

To simplify the considerations, we shall neglect nonlocal terms. In our approch, this 
is justifiable since if we are able to solve the singularity problem in a purely local theory, 
we expect that the nonlocal terms (which are inevitable for example because of particle 
production) will not drastically change the asymptotic behavior of our theory because of its 
special properties (see Section 5). 
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Key to the analysis is the assumption about validity of the background field approxima- 
tion for the gravitational field up to high curvatures. Such an approximation will only be 
justified if the quantum fluctuations around this metric are sufficiently small. If the gravi- 
tational field is asymptotically free at high curvatures (see Section 3), we can hope that this 
approach will be valid. As we shall see, there are features in our theory which indicate that 
this will really be the case. 

For the moment we shall ignore matter (later we will show that the presence of matter 
does not change the solutions at high curvatures in an important way). Thus, our starting 
point is the effective action 

s g = Jf(r,r^r^,c 2 ,...) V=gd 4 x. (2.2) 

The usual Einstein theory in the absence of matter has only one solution - Minkowski 
space - for a homogeneous and isotropic Universe. Any non-Einstein theory of gravity 
gives rise to fourth (or higher) order equations of motion and hence to a large number of 
cosmological solutions. In general, the singularity problems of such a theory are much worse 
than in Einstein gravity. A simple example is R 2 gravity 

F(R) = R + aR 2 (2.3) 

which is conformally equivalent 19 ) to Einstein gravity plus scalar field matter and which 
hence has many isotropic singular solutions (even without matter). Thus, the theory we are 
looking for must be a very special higher derivative gravity model. 

We wish to construct an effective action for gravity in which all homogeneous and 
isotropic solutions are nonsingular and in which all curvature invariants are limited (in 
Section 5 we will indicate how to extend our analysis to anisotropic models 13 )). To motivate 
our construction, it is useful to keep in mind ways of writing the action for two well-known 
physical theories in which certain physical quantities are bounded: special relativity and the 
Born-Infeld theory of electromagnetism 20 ) . 

To impose bounds on physical quantities in an explicit manner, it is convenient to employ 
a Lagrange multiplier technique proposed by Altshuler 12 ) . To explain how this technique 
works, we first consider the simple example of point particle motion. We start with the 
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action for a nonrelativistic particle of mass m and world line x(t). We demonstrate how to 
explicitly implement the limitation on the particle velocity, and in particular how to obtain 
the action for point particle motion in special relativity. The nonrelativistic action with 
which we start is 



'old 



= m J dt^x 2 . (2.4) 



In order to construct a new theory with bounded velocity we introduce a "Lagrange mul- 
tiplier field" (pit), which couples to some function of the quantity whose value we want to 
limit, and a potential V{4>) for this field: 



m I dt 



\x 2 + (px 2 -v(</>) 



(2.5) 



Let us stress that <fi is not a dynamical field. Provided that dV/d<f> is bounded, the constraint 
equation (i.e. the variational equation with respect to 0) ensures that x is bounded. In order 
to obtain the correct Newtonian limit for small x and small , V(<f>) must be proportional to 
4> 2 as |0| — > 0. One of the simplest potentials which satisfies the above asymptotic conditions, 

^> = TW < 2 - 6 > 

leads to special relativity. In fact, eliminating the Lagrange multiplier using the constraint 
equation and substituting the result into (2.5) yields (up to a constant term which does not 
effect the equations of motion) the relativistic point particle action 



m 



J dt v 7 ! - x 2 . (2.7) 



Let us return to the theory of gravitation. In the notation of the above example, the 
"old" theory will be given by the Einstein action. In order to implement the LCH, we wish 
to impose restrictions on some curvature invariants I\ , I2 , ■ ■ ■ , I n m an explicit manner. 
The general form of a higher derivative local modification of the Einstein action involving 
the invariants I\ , . . . , I n is 

S 9 = --[^q J [R + F(h,I 2 ,...,In)} V=gd A x, (2.8) 
where F is some function of the invariants I\ , ... , I n . 
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By introducing Lagrange multiplier fields (j>i(t) , ...,(f> n (t), the above action can be 
rewritten as 

Sg = ~ i6^G / \. R + ^M^) + ••• + MM + V {</>!... (M] V=gd 4 x, (2.9) 

where fi{U) are functions we can choose as we want. The actions (2.8) and (2.9) are equiv- 
alent provided that the potential V{(f>\ , . . . , <p n ) satisfies the following partial differential 
equation 

dV 



-E*^+^-*.) = *"(a-(S.-^( 

This follows immediately by using the constraint equations for (2.9) 



d(f) r , 



(2.10) 



dV 

f^k)^— i = l,...,n. (2.11) 

We see from the constraint equations (2.11) that by appropriate choice of the functions 
fi and V we can implement bounds on the invariants I\ , . . . , I n . Variation of the action 
(2.9) with respect to g^ v yields the field equations.. 

First, we try to construct the simplest theory in which the LCH is realized. At least for 
simple models (like the isotropic Universe) it is natural to choose as one of the invariants 

h = R- \/?>{±R^R^ - R 2 ) l/2 (2.12) 

since for a homogeneous, spatially flat Universe it is equal to 12H 2 . This invariant will be 
used to limit the curvature by some (e.g. Planckian) value. The second invariant I2 we will 
take on such a form as to implement in the theory the condition that in the asymptotic 
regions all of the solutions evolve to de Sitter. The simplest way to do this is to pick I2 such 
that I2 = only for de Sitter space (Minkowski space is included as a special case), and to 
make sure that 

I 2 ^0 as |0 2 |^oo. (2.13) 

For homogeneous and isotropic space-times, it can be shown that 

h = AR^R» v -R 2 (2.14) 
is a good choice, since h = only for de Sitter space. Note that in general, I2 is positive 



9 



semidefinite. However, for inhomogeneous and anisotropic space-times (e.g. when C 2 ^ 0), 
the above form of I2 is insufficient to single out de Sitter space as an asymptotic solution. 
This is obvious from considering the Schwarzschild metric for which J 2 = 0. Hence, in the 
general case we 13 -* should add to (2.14) terms which depend on C 2 and vanish for conformally 
flat space-times. 

However, for a homogeneous and isotropic Universe it is (as we will show) sufficient to 
consider the action in the following general form: 

Sg = ~ltc J [ jR + ^^( / i) + ^/2(/ 2 ) + U(0 1 ,0 2 )] ^d 4 x. (2.15) 
The variational field equations which follow from (2.15) are 

Rs-bsR-bsV=(^ dflY ^ dflY 



V ^2 J 3 
,0 / of \ ,a 



+ n ( 02 ||«) • _ 2 a Ji R J 802 a M ^ 



g /2 „ na , 1 



+ 2 2 ^ J R J R« + -^0 2 / 2 , 



and the constraint equations are 

=- 



(2.17) 



We will simplify the theory further by assuming a factorizable potential 



^(0i,0 2 ) = U 1 (0 1 )+V2(0 2 ). (2.18) 



The asymptotic conditions on the potentials V\ and V 2 follow from demanding that the 
theory reduces to the Einstein theory at small curvatures, and that the LCH is realized. The 
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first condition yields 



V r i (0i)~0? |&|«l;i = l,2 



(2.19) 



In order to limit R explicitly we can try a potential which to leading order takes the form 



Vi(0i 



|0l| » 1, 



(2.20) 



and to obtain de Sitter solutions in the asymptotic regions we need a potential which at 
large 02 increases less quickly than 02- We assume an asymptotic form: 



V2(<f>2) ~ const 1 02 1 > 1 ■ 



(2.21) 



In this case, provided ^2(^2) - ► as I2 — > 0, the constraint equation (2.17) implies that 
J2 — > as |02 1 — > 00, and we have a chance of realizing the LCH, provided that the evolution 
of the scalar fields 0i and 02 is appropriate, a question which needs detailed investigation. 

To conclude this section we will write down equations (2.16) and (2.17) explicitly for 
a homogeneous and isotropic metric with scale factor ait) in the contracting phase (i.e. 
H < 0): 



ds 2 = dt 2 - a 2 (t) 



1 



dr z + r 2 d6 2 + r 2 sin 2 9 d<p 



1 — kr 2 

We choose simple functions f\ and fx 



(2.22) 



h{h) = h 
f 2 (h) = -Vh 



(2.23) 



Thus, our final action takes the form 



(2.24) 



As is well known from the derivation of the Friedmann-Robertson- Walker equations in 
Einstein gravity, the only independent equation of motion is the — equation. In our 
case, we have in addition the constraint equations (2.17). The full set of equations can be 
obtained by inserting the metric (2.22) into (2.16) and (2.17). 
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The resulting 0i , 02 and — equations are 



k 1 



H-±=-j=VZ (2.26) 



- \{Vi + U 2 ) + 3# 2 (1 - 20i) + 3 ^(40i + 1) = V3H (^0 2 + 3#0 2 - ^<hj (2-27) 

Another way to obtain the same equations is to substitute the ansatz (2.22) with goo = N(t) 2 
into the action (2.24) and to vary with respect to N ,<p\ and 02 (see e.g. Ref. 21). Adding 
to the system matter with action 

S m = J L m ^d A x (2.28) 
where L m is the matter Lagrangian, only leads to an additional term 

Y°Pm (2-29) 

on the left hand side of the — equation. 

In the following sections, we shall show that all solutions of the above equations are free 
of singularities. 



3. Nonsingular Universe without Limiting Curvature 

Since our goal is primarily to construct a nonsingular Universe model and only secon- 
darily to limit the curvature, we first consider a simple model in which the 0i field is absent. 
In this case, it is easier to discuss our techniques of analysis. 

We will show that for this model all solutions for a collapsing Universe are nonsingular 
and asymptotically approach de Sitter solutions. However, there is no general (i.e. solution 
independent) bound on the effective cosmological constant of the de Sitter period. 
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In this section we set 02 = an d Vi = V. The equations of motion are given by (2.26) 
and (2.27). Let us first consider a spatially flat (k = 0) collapsing model without matter. In 
this case, the equations of motion are 



= - 3H(f) + V3H ^-V . (3.2) 

2v 3 H 

The phase space of this model is the two dimensional (0 , H) plane. The phase space 
trajectories can be understood by considering dH/d<f> (determined immediately from (3.1) 
and (3.2)) 

iH V ' {-ZH, + ^H-^v)-\ (3.3) 



d<p y/Y2 

From (3.3), it follows that provided that V(<f>) is bounded at large (as postulated in 
(2.21)), then as tends to infinity, H approaches a finite value, i.e. for any solution, the 
effective cosmological constant in the large region is bounded. In this case, it follows from 
(3.2) that in a collapsing Universe, for large 

0(t)~e 3 ^l t . (3.4) 

Our choice of invariant I2 has led to the conclusion that the asymptotic de Sitter solutions are 
attractor solutions. This conclusion holds independent of the specific choice of the potential 
V(0), as long as the asymptotic condition (2.21) is satisfied. 

From (3.4) it follows that all solutions for a contracting Universe are free of singularities. 
It takes infinite time to reach = oo. 

To concretize the consideration, we consider a simple potential which satisfies the asymp- 
totic conditions (2.19) and (2.21): 

V = y/n #0^2, (3-5) 

where Hq is a constant (in the model with limiting curvature discussed in Section 4, Hq sets 
the scale of this limiting curvature). 
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The phase space trajectories (0(i) ,H(t)) in a collapsing Universe are shown in Fig. 2. 
The numerical results were obtained using the specific potential (3.5). However, as discussed 
above, the main features of the diagram depend only on the asymptotic properties. 

First we note that there is only one singular point (0 = H = 0) in the phase plane. This 
point is 

(0, if) = (0,0) (3.6) 

and corresponds to Minkowski space-time. 

There are two classes of trajectories which are asymptotically de Sitter. Those starting 
at large positive values of go off to = oo, reaching their asymptotic value of H from 
above (i.e. H < 0). Those starting with large negative values of tend to = — oo with 
H > 0. 

For small values of H and we can use the asymptotic condition (2.19) on V(4>) to 
conclude that there are periodic solutions about Minkowski space. In this limit, the basic 
equations (3.1) and (3.2) become 

(3 - 7) 

j4^ gE o^^, (3.8) 
V y/3 H H/H ' 1 J 

where for V(<f>) we have inserted the general asymptotic form 

V{<j>) ~2v / 3/f o 2 2 , (3.9) 

valid for small 0. The numerical factor 2\/3 has been inserted to eliminate numerical con- 
stants in the following equations. 

It is convenient to introduce a rescaled time 

r = H t (3.10) 

and a dimensionless measure of H: 

y = H/H (3.11) 
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With d/dr denoted by a prime, equations (3.7) and (3.8) become 



V = -20 

, VSy 2 -^ 2 (3-12) 



y 



To see the oscillatory nature of the solutions, we introduce radial and angular coordinates r 
and ip 

<p = r sin if) 

y = — 3 1 r{l — cosip) 
The resulting equations for r and ip are 

(3.14) 

r' = 0. 

where the frequency is cj = 2 • 3 1 / 4 . The corresponding solutions oscillate with frequency 
given by Hq (which we expect to be Planck scale) about Minkowski space. 

Based on the preceeding discussion of asymptotic solutions we see that there is a separatrix 22 )| 
in phase space dividing solutions which tend to <fi = oo from those which oscillate or tend to 
= — oo. We observe that for large \H\, the separatrix will asymptotically (and from the 
right hand side on Fig. 2) approach the line of turning points given by d<p/dH = 0. From 
(3.1), it follows that for large \H\ the turning points lie at 

(3.15) 

For small values of and H, the separatrix is well to the right of the line of turning points 
given by 

0~3 1 / 4 ^. (3.16) 

The above analysis of the phase space trajectories is an indication that in our theory, 
Minkowski space is stable towards homogeneous perturbations. As long as the initial values 
of | H | , |0| and <f>/\H\ are small, a solution starting close to Minkowski space will remain close 
for all times. The issue of stability of Minkowski space towards inhomogeneous perturbations 
is an important unsolved problem. 
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We stress again that all the general features of the phase space analysis are true for any 
potential V{4>) which satisfies the required asymptotic conditions (2.19) and (2.21). However, 
the results depend crucially on the choice of the invariant Ii- 

Next, we include hydrodynamical matter with energy density 

p m (t) = ca(t)- n , (3.17) 

where n = 3 for dust and n = 4 for radiation. For the moment, we keep to a collapsing 
spatially flat model. In this case, equation (3.1) is unchanged while equation (3.2) becomes 

<j> = -3H<f> + y/3H t^V- ^^ca(t)- n . (3.18) 

With matter, phase space is three dimensional - the third dimension being a(t). In 
Fig. 3, we show the projection of some of the trajectories onto the (<j>(t) , H(t)) plane, for 
potential V(<f>) given by (3.5). All trajectories have 8nG c = 1 and a (to) — 10, to being the 
initial time. The main impression is that the trajectories look very similar to those without 
matter in the asymptotic region. We shall now explain why this is the case. 

First, we note that as \<f>\ — > oo, the solutions approach de Sitter space since H — > 0. 
Hence, 

a(t) ~ e-WQ-* ) a(t Q ) . (3.19) 
Next, we combine (3.1) and (3.18) to obtain for |0| >> 1 

dH V / , 8nGc 



_ . :;// n(t)-A . (3.20) 

Our model incorporates a very important feature: in the asymptotic de Sitter region matter 
does not have an important effect on the geometry. The effective gravitational constant which 
describes the influence of matter on the geometry goes to zero as space-time approaches de 
Sitter space. In this sense the model is asymptotically free. 
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Some understanding of asymptotic freedom can be obtained by solving the and H 
equations of motion (3.1) and (3.18) in the asymptotic region |0| >> I. Equation (3.18) 
becomes 

~ 3|#|0 + -^-a(t )- n e nmt - to) (3.21) 
\H | 

(where we have incorporated the factor 8ttG/\/3 into the definition of c). From (3.21) it 
follows that <p(t) is a linear combination of the homogeneous solution (3.4) and (assuming 
that H ~ const ) the inhomogeneous contribition <pi(t) 

Mt) = z^m a ^r U (e^K'-'°) - l) (3.22) 



n 



For dust {n = 3), both the homogeneous and inhomogeneous terms grow at the same rate, 
and the coefficient of the inhomogeneous term is smaller. Hence, matter does not effect even 
the time dependence of the phase space trajectories. For radiation (n = 4), 4>i(t) grows faster 
than (3.4). At sufficiently late times, therefore, it will dominate. In this period, however, we 
can (for potential (3.5)) solve the H equation (3.6) to obtain 

3 



H(t) ~ H(t!) - [ ^MEL \ a(t f n e ~ Sn \ H ^-^ , (3.23) 

3n\H\ \ c / 

(where t\ is some time 3> to we U m to the asymptotic region) which shows that the presence 
of matter does not effect the final value of the curvature when starting the evolution in the 
asymptotic region. 

For small |0|, the presence of matter does have a significant effect on the phase space 
trajectories. As a(to) decreases (or, equivalently, c and thus the matter energy density 
increase), the distortions of the trajectories increase, as can be seen by comparing Figs. 3 
and 4. Figure 4 corresponds to a matter energy density which is ten times larger. 

Finally, we consider the effects of spatial curvature. In this case, equations (3.1) and 
(3.2) generalize to (see (2.26) and (2.27)) 

"-^""4 (3 - 24) 



where the constant c is as in Eq. (3.21). 
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In the case of the potential (3.5) and for c = 0, some resulting phase space trajectories 
projected onto the ((f)/ H) plane are shown in Figs. 5 and 6. For the trajectories of Fig. 6, 
the initial value of a(t) was chosen to be ten times smaller than in Fig. 5. Hence, the effects 
of curvature are more pronounced. 

Consider a sample trajectory of Fig. 5. It starts out with large initial value of a. The 
trajectory tends towards \<p\ » 1 and H — > 0, as in the case k = 0. Since a{t) is now 
decreasing almost exponentially, the role of curvature increases. At a critical value of <j>, the 
value of H becomes 0. This will occur when 

Hence, the smaller the initial value of a(t), the earlier (3.25) will be satisfied (compare Figs. 
5 and 6). At a similar time, the curvature terms also start to dominate in Eq. (3.24). 
Therefore, as is obvious from the k dependent terms in (3.24), <fi(t) will rapidly decrease, as 
will \H(t)\. At some finite and negative value of 0, H(t) vanishes. Thereafter, the Universe 
reexpands. The evolution of this model for small a(t) resembles a de Sitter bounce. 

Note that all solutions are nonsingular. In particular, the solutions can be integrated 
through the point when H = (when terms on the right hand side of (3.29) become infinite). 

In conclusion, we have constructed a higher derivative modification of Einstein's theory 
in which all homogeneous and isotropic solutions are nonsingular. Without curvature (i.e. 
for k — 0), the solutions either are periodic about Minkowski space or else converge to a 
k = de Sitter solution. For k ^ the solutions which do not remain close to Minkowski 
space go through a de Sitter bounce and are future extendible to t = oo. In addition, we 
have shown that our model is asymptotically free in the sense that the effective coupling of 
matter to gravity goes to zero as the curvature increases. 
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4. Nonsingular Universe with Limiting Curvature 

Now we turn to the discussion of the full model in which the LCH is implemented, 
the model given by the action (2.24), in which for a homogeneous and isotropic metric the 
equations of motion reduce to (2.25 - 2.27). We include hydrodynamical matter with energy 
density given by (3.17). 

In the general case (k ^ and c ^ 0), the phase space of the model is three dimensional: 
0i(£) ) 02 (t) and a(t). For k — and c = 0, the dependence on a(t) drops out and the phase 
space can be reduced to the two dimensional 0i/02 diagram. The first order equations 
of motion in phase space are found by combining equations (2.25 - 2.27). To derive the 
equation for <f>i(t), we differentiate (2.25) with respect to t and use (2.26) to substitute for 
H to obtain 

0i = -4V^. (4.1) 

The equation of motion for 02 is (2.27): 

2 = -ZH<h + ^02 + (SH 2 (1 - 20! ) + 3^(40! - 1) - i(Vi + V 2 ) - (4.2) 

where H can be expressed in terms of 0i and a via (2.25). From (4.1), (4.2) and (2.25) it is 
obvious that for k = c = the a(t) dependence disappears. 

In the case k = c = we may use (2.25) to get 



d02 V{' 



-v / 30 2 + (l-20i)-| 7 (Vi + V2) 

v i 



(4.3) 



#1 4V^ 

the key equation for the following phase space analysis. 

For all potentials Vi(0i) and 1^2 (02) with asymptotical behavior 

Vicxtf 0i«l (4.4) 

and 

Vi oc 0i - ln0i + 0(— ) 0i > 1 (4.5) 
01 

V 2 oc const + O(-J-) 02 > 1 (4.6) 

02 

the phase diagrams have the same features as depicted schematically in Fig. 8 for spatially 



19 



collapsing Universes without matter. The numerical solutions depicted in Fig. 7 were ob- 
tained for the particular choice of potentials (which satisfy the asymptotic conditions of (4.4 



The presence of the logarithmic term in (4.4) will be justified shortly. 

We can identify four classes of trajectories. Note that by (2.26), |02| - ► oo implies that 
the evolution approaches de Sitter space. The first class of trajectories starts in the de 
Sitter phase at 02 —> — oo and evolves to de Sitter at 02 — > oo. For small initial values of 
01, trajectories starting at 02 = — oo reach a turning point and return to 02 = — oo. The 
third class of trajectories are periodic solutions about Minkowski space-time (0i = 02 = 0). 
Finally, trajectories starting with small 0i and 0i/02 with 02 positive evolve towards de 
Sitter solutions at 02 = oo. There are two separatrices dividing phase space into regions 
corresponding to the four above classes (see Fig. 8). 

Note that in order to prevent solutions starting with 0i >> 1 and 02 ~ from escaping 
to 0i = oo at 02 < 1 in finite time - such solutions would violate the LCH and would lead 
to singularities in higher order curvature invariants - it was necessary to add the logarithmic 
correction term to Vi(0i). 

Phase space is the half plane 0i > 0. Negative values of 0i are unphysical since by (2.25), 
and using the small 0i asymptotic form of Vi(<f>i), they would correspond to imaginary 
values for H(t). This half plane can be divided into four regions: in Region A, 0i >> 1 and 
|02| >> 1, in Region B, 0i >> 1 and |02| << 1, in Region C, 0i << 1 and |02| << 1 and 
in Region D, 0i << 1 and |02| >> 1. We will analyze the phase space trajectories in each 
of the above regions, focusing on three features: the asymptotic expressions for <i02/<i0i 
(which give the tangent vectors to the trajectories), the separatrices, and the equations 
for the trajectories. To concretize the discussion, we use the potentials (4.7) and (4.8). 
However, except in Region B, the asymptotical solutions are independent of the specific 
choice of potentials. 



- 4.6)) 




(4.7) 



U 2 (0 2 ) = 2y/3H$ 



1 + 02 



(4.8) 
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In Region A, Equation (4.3) becomes 

d(f)2 V3 0| 



02 



(4.9) 



#1 ~ 4 02 

The direction of the tangent vectors is sketched in Fig. 8. Arrows indicate the direction of 
increasing time and are obtained by inspecting (4.1) and (4.2) directly. By inspecting the 
tangent vectors it is clear that all solutions in the upper region 02 > 1 quickly approach de 
Sitter space (|02| — > oo implies de Sitter space). In the lower region 02 < —1 there are two 
domains separated by a separatrix which for 0i ^> 1 and |02| 1 is close to the line of 
turning points where dfo/dcj)! = 0, its equation being given by 

(see Fig. 8). To the right of the separatrix, trajectories correspond to solutions starting 
out in a de Sitter phase. To the left of the line given by (4.10), we have dfo/dcfii > and 
trajectories go off to de Sitter space at 02 — > — oo. In all cases, de Sitter space is reached 
at finite 0i values. This is seen by explicitly integrating (4.9). In the region where the first 
term on the right hand side of (4.9) dominates we have 

1 ~c-^0 2 - 3 , (4.11) 

while in the domain where the second term dominates the approximate solution is 

01 * 9-^— (4.12) 

(c is a constant of integration). Note that all of the solutions starting in Region A start in 
de Sitter space and end up in de Sitter space. 

In Region B, the tangents in phase space are given by 



which integrates to 

0i = cexp|^i=02| (414) 

The tangent vectors are again sketched in Fig. 8. From (4.14) it follows that trajectories 
leave Region B at a finite value of 0i. They enter Region A and hence asymptotically 
approach de Sitter space. 
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In Region C, Equation (4.3) becomes 

d(f)2 a/3 



1 - 30i - 



1 



4£ 

12 0i 



(4.15) 



The separatrix in the upper half planes is close to the line of turning points dfo/dcpi = for 
large values of 4>f 



02 ^ ±12 1/4 0l /2 (l-30; 

Where the first term dominates, the trajectories obey 

01 ^ —j= (Pl + C 



,1/2 



(4.16) 



(4.17) 



From the sketch of Fig. 8 it is clear that the trajectories which pass through 0i = 02 = 
with 02/01 (01 = 02 — 0) not too large correspond to periodic motion about Minkowski 
space. This - as in the model of Section 3 - is an indication that Minkowski space is stable 
in our theory towards homogeneous perturbations. 

Finally, in Region D the equation for the tangent vector is 



d<pi 



>/3, 3 



a/302 + \ 

2a/30i 



(4.18) 



There is a separatrix which is (for large (pi) approximately described by 

1 



02 



60i 



(4.19) 



To the right of this line, the trajectories are given by 

2 , 
01 - c - - 02 , 



to the left by 



(4.20) 



(4.21) 



In conclusion, all solutions are either periodic about Minkowski space or are asymp- 
totically de Sitter. All solutions can be extended to t = ± oo, and hence there are no 
singularities. 
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As in the model of Section 3, including matter does not effect the asymptotic solutions. 
The coupling between matter and gravity is asymptotically free also in the theory with action 
(2.24). However, including matter changes the nature of solutions starting near Minkowski 
space. These solutions now approach de Sitter space (see Fig. 9). This result is not surpris- 
ing, since also in Einstein gravity, Minkowski space is not a solution of the field equations in 
the presence of matter. 

The projection of some phase space trajectories onto the (0i , 02) plane in a model with 
k 7^ but c = are shown in Fig. 10. Like in the single field model of Section 3, the 
trajectories initially evolve as for k = towards de Sitter space. Hence, for finite 0i , 02 
becomes very large. Eventually, however, the curvature terms become important, 02 reaches 
a turning point and rapidly (within time period Hq 1 ) relaxes to zero (for finite value of 0i). 
As is obvious from Fig. 11, the rapid decrease in 02 corresponds to the de Sitter bounce 
during which H changes sign. 

5. Conclusions and Discussion 

We have constructed a theory of gravity in which all homogeneous and isotropic solutions 
(not only special solutions as in some other mo dels 23 )) are nonsingular, regardless of the 
matter content of the Universe. Our effective action for gravity contains higher derivative 
terms which modify the Einstein action at high curvatures. Such terms are expected to be 
important near the Planck curvature in any fundamental theory such as quantum gravity or 
string theory. 

Most higher derivative gravity theories have much worse singularity properties than 
Einstein gravity. We use a particular construction based on implementing the "Limiting 
Curvature Hypothesis" to obtain a class of models without singularities. We discussed two 
models, one in which all curvature invariants are bounded and all solutions except those 
periodic about Minkowski space asymptotically approach de Sitter space (Section 4), and a 
simpler model without limiting curvature (Section 3). 

The theory presented in this paper is "asymptotically free" in the sense that the coupling 
of matter to gravity goes to zero as the curvature approaches its limiting value (similar 
features have been discussed by Linde 24 ) under the name "gravitational confinement"). 
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When applied to an expanding Universe, our theory implies that it started out in a de 
Sitter phase with scale factor a(t) = e Ht (for k = 0) or else (for k — 1) it emerged from a 
de Sitter bounce. In particular, there was a period of inflation driven by gravity This is no 
surprise as it is well known 25 ) that higher derivative gravity theories often produce inflation. 

Note that the property of asymptotic freedom might also justify using the effective action 
approach to gravity until Planck curvatures. Asymptotic freedom will also play an important 
role in controlling nonlocal terms. For example, nonlocal terms due to particle production 
may be expected to vanish in the asymptotic regions of phase space. 

Our action is constructed by adding two Lagrange multiplier terms (and their correspond- 
ing potentials) to the Einstein action. Each Lagrange multiplier is coupled to a curvature 
invariant. The role of the first Lagrange multiplier is to limit the curvature, the role of the 
second one (02) is to force space-time to be de Sitter at large curvature. For a homogeneous 
and isotropic model, it was sufficient to couple 02 to the invariant I2 = AR IJLV R tJ,v — R 2 , since 
in this case h = singles out de Sitter space. 

However, for an anisotropic cosmology, we must extend the invariant I2 by including 
a term which effects the anisotropy. In a subsequent paper 13 ) (see also Ref. 26) we show 
that h = +R„uR^ -R 2 + C 2 is an appropriate invariant. This invariant also works for a 
spherically symmetric metric. Thus, in a model like the one presented here, but with the new 
I2, we are able to show that there will be no singularities inside the black hole horizon 13 ). 

Open questions include the generalization of our model to general inhomogeneous met- 
rics, and a full stability analysis. 
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Figure Captions 

Fig.l : Penrose diagram of an eternal black hole in Einstein gravity (a) and in the nonsingular 
Universe theory (b). The singularities (S) are replaced by de Sitter phases (dS) which 
couple to Friedmann Universes (FRW). The horizons (H) are not affected. 

Fig.2 : Phase-space diagram (0, H) - arrows indicating the direction of time evolution - for 
a spatially flat Universe without limiting curvature and with no matter (k — c — 0). 
Generated using the potential (3.5). 

Fig.?> : A projection onto the (0, H) plane of the three dimensional phase-space diagram 
((f), H, a) - arrows again indicating direction of time evolution - for a spatially flat 
Universe without limiting curvature but with matter (k = 0,c ^ 0). Generated using 
the potential (3.5) with initial condition a (to) = 10. 

Fig A : Phase-space diagram as in Fig. 3, but with a (to) = 1- Therefore, the initial matter 
energy density is larger than for the trajectories of Fig. 3. 

Fig.5 : A projection onto the (0, H) plane of the three dimensional phase-space diagram 
((p, H, a) in a closed (k — 1) Universe without limiting curvature and in the absence 
of matter (c = 0). Potential (3.5) was used, and a(to) = 10 was chosen as initial 
condition. 

Fig.6 : Same as in Fig. 5, but with initial condition a(to) = 1. Notice the different scales on 
the axes. 

Fig.l : Phase-space diagram for the spatially flat (k = 0) Universe with limiting curvature 
based on the potentials (4.7) and (4.8). There is no matter (c = 0). 

Fig.8 : Sketch of the generic phase-space diagram for a two field model with k = c = 
and potentials satisfying the asymptotic conditions (4.4) - (4.6). Lines with arrows 
indicate phase space trajectories, arrows pointing in the direction of increasing time. 
Separatrices are shown as dashed lines. With A, B, C and D we denote the asymptotic 
regions of phase space discussed in the text. 

Fig. 9 : A projection onto the (0i,02) plane of the three dimensional phase-space diagram 
(01, 02, a), for a two field model which is spatially flat but contains matter. The 
potentials used are (4.7) and (4.8). 

Fig. 10 : The same for a two field model without matter but including spatial curvature (k ^ 0). 



27 



Fig.ll : Trajectories in the (h, H) plane for the same model as in Fig. 



28 



